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Abstract — Percolation theory has become a useful tool for 
the analysis of large-scale wireless networks. We investigate 
the fundamental problem of characterizing the critical density 
\^c^ for c{-dimensional Poisson random geometric graphs in 
continuum percolation theory. By using a probabilistic analysis 
which incorporates the clustering effect in random geometric 
graphs, we develop a new class of analytical lower bounds for the 
critical density \\f^ in c{-dimensional Poisson random geometric 
graphs. The lower bounds are the tightest known to date. In 
particular, for the two-dimensional case, the analytical lower 

("21 

bound is improved to Ac > 0.7698.... For the three-dimensional 
case, we obtain A^^' > 0.4494.... 

I. Introduction 

Recently, percolation theory has become a useful tool for 
the analysis of large-scale wireless networks [l]-[4]. A per- 
colation process resides in a random graph structure, where 
nodes or links are randomly designated as either "occupied" or 
"unoccupied." When the graph structure resides in continuous 
space, the resulting model is described by continuum percola- 
tion [5]-[10]. A major focus of continuum percolation theory 
is the d-dimensional random geometric graph induced by a 
Poisson point process with constant density A. A fundamental 
result for continuum percolation concerns a phase transition 
effect whereby the macroscopic behavior of the system is very 
different for densities below and above some critical value 
Ac^''. For A < Ac'''' (subcritical), the component containing 
the origin contains a finite number of points almost surely. 
For A > Ac''' (supercritical), the component containing the 
origin contains an infinite number of points with a positive 
probability [7]-[10]. 

Naturally, the characterization of the critical density Ac''' 
is a central problem in continuum percolation theory. Un- 
fortunately, the exact value of Ac''' is very difficult to find. 
For two-dimensional random geometric graphs, simulation 



cases. 

In this paper, we give a new mathematical characterization 
of the critical density aI''' for Poisson random geometric 
graphs in d-dimensional Euclidean space, where d > 2. 
We develop an analytical technique based on probabilistic 
methods [13] and the clustering effect in random geometric 
graphs. This analysis yields a new class of lower bounds: 
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where V^'^'' is the volume of a d-dimensional unit sphere 
and Ct''' is the t-th order cluster coefficient it > 3) for d- 
dimensional Poisson random geometric graphs, which we will 
define later. This class of analytical lower bounds are the 
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tightest known to date. In particular, by evaluating C3 in 

closed-form, the analytical lower bound for two-dimensional 

(2) 

Poisson random geometric graphs is improved to Ac > 

0.7698.... For three-dimensional Poisson random geometric 

(3) 

graphs, we obtain the analytical lower bound Ac > 0.4494.... 
By successively evaluating Cj''' for t > 4, we can obtain 
tighter lower bounds on Ac ' . 

II. Random Geometric Graphs 

In wireless networks, a communication link exists between 
two nodes if the distance between them is sufficiently small, so 
that the received power is large enough for successful decod- 
ing. A mathematical model for this scenario is as follows. Let 
II • II be the Euclidean norm, and / be some probability density 
function (p.d.f.) on M''. Let Xi,X2,...,X„ be independent 
and identically distributed (i.i.d.) d-dimensional random vari- 
ables with common density /, where denotes the random 
location of node i in W^. The ensemble of all the graphs 
with undirected links connecting all those pairs {xi,Xj} with 



studies show that Ac ~ 1.44 [11], while the best analytical ||xi — Xj || < r, r > 0, is called a random geometric graph [8]. 

In the following, we focus on random geometric graphs 
with Xi,X2, ...,X„ distributed i.i.d. according to a uniform 
distribution over a given d-dimensional box A = [0, y/nJXf' 



bounds obtained thus far are 0.696 < A^^' < 3.372 [6], 
[7]. Recently, in [12], the authors reduce the problem of 

(2) 

characterizing Ac to evaluating numerical integrals using a 



mapping between continuum percolation and dependent bond 

percolation on lattices. By Monte Carlo methods, they obtain 

(2) 

numerical bounds 1.435 < Ac < 1.437 with confidence 
99.99%. Unfortunately, the bounds obtained in [12] are not in 
closed-form and cannot be generalized to higher dimensional 



with A > 0. We denote such graphs by G{X^f^ ; r). 

Let A = \A\ be the d-dimensional Lebesgue measure (or 
volume) of A. An event is said to be asymptotic almost sure 
(abbreviated a.a.s.) if it occurs with a probability converging 
to 1 as n ^ 00. Consider a graph G = (V, E), where V 



and E denote the set of nodes and links respectively. Given 
u,v ^ V, we say u and v are adjacent if there exists an link 
between u and v, i.e., (u, v) € E. In this case, we also say 
that u and v are neighbors. 

A. Preliminaries 

In G(A'i'*''; r), let the location Xj of node j be given. A 
second node i is randomly placed in A according to the 
uniform distribution /. There exists a link between these two 
nodes if and only if node i lies within a sphere of radius r 
around Xj. Let this spherical region be denoted as ^(xj ), then 
the probability for the existence of a link between i and j is 
given by 

P/mfe(x,) = / /(y)dy. (2) 

Since the underlying distribution is uniform, the probability 
Piinki^j) depends only on the volume of the intersection 
between A and the node coverage volume ^(x^). Throughout 
this paper, we ignore border effect^ As a consequence, 
Piinki^j) is independent of Xj, and thus independent among 
all the links: 

Plink — -: — , (3) 

y\ Ti 

where V'^'^^ is the volume of a d-dimensional unit sphere 
V^'^'> = Y^TT^ and r(-) is the Gamma function r(a;) = 
t^-^e-'dt. 

It follows that in G(A'^'*^; r), the probability that the given 
node j has degree A:,0<fc<n — 1, is given by the binomial 
distribution: 

Pk^{^l^)pLk[^-PunkT-''^. (4) 

Thus, the mean degree for each node is 

^x = E[k]^{n~l)Punk^^- '— . (5) 

n 

Note that as n and A both become large but the ratio n/A = 
A is kept constant, each node has an approximately Poisson 
degree distribution [8], [14] with an expected degree 

lim . Ay(^V. (6) 

n — *oo 77, 

As 77 — > cx) and A oo with 77/^4 = A fixed, G{x!!i'^\r) 
converges in distribution to an (infinite) random geometric 
graph G{Ti}-^^\r) induced by a homogeneous Poisson point 
process with density A > 0. For such graphs, we have the 
following lemma. 

Lemma 1: Suppose G{H^^^]r) is a random geometric 
graph induced by a homogeneous Poisson point process with 
density A > 0. Let A' be any subset of M'* with \A'\ < oo, 

'More rigorously, we may use a torus instead of a box for A. Asymp- 
totically, as n — » oo and A — > oo with n/A = X fixed, uniform random 
geometric graphs on the torus are the same as those in the box. 




Fig. 1. Calculation of cluster coefficient C'^' 

then the subgraph G' contained in A' has a finite number of 
nodes a.a.s. 

Proof: The proof is straightforward and omitted here. □ 

Due to the scaling property of random geometric graphs [7], 
[8], in the following, we focus on G{H^f^; 1) and Gixi'^^; 1). 

B. Cluster Coefficients 

An important characteristic of random geometric graphs 
is the clustering effect. Here, if node i is close to node j, 
and node j is close to node k, then i is typically also close 
to k. In the following, we use the cluster coefficients to 
precisely characterize the clustering property. This turns out 
to be the key to deriving new bounds for the critical density 
in continuum percolation. 

Definition 1: Given distinct nodes k in G'(A'i'''; 1), the 
cluster coefficient C'^'^^ is the conditional probability that nodes 
i and j are adjacent given that i and j are both adjacent to 
node k. 

The calculation of C^"^^ for two-dimensional random geo- 
metric graphs is illustrated in Figure 1. To determine C^^^, as- 
sume both nodes i and j lie within A{'x.k), then the conditional 
probability that nodes i and j are also adjacent is equal to the 
probability that two randomly chosen points in a circle with 
radius 1 is at most distance 1 apart. In other words, given the 
coordinates of x^ and x^, the probability that there is an link 
between i and j is equal to the fraction of yl(xi) that intersects 
A{yik)- By averaging Xj over all points in ^(x^), the cluster 
coefficient can be found as C^^^ = 1 - « 0.5865... [15]. 

The cluster coefficient C*^''^ reflects the triangle effect in 
random geometric graphs. To further capture the cluster effect, 
we now generalize the notion of cluster coefficients for more 
than three nodes. 

Definition 2: For t > 3, suppose vi, . . . , vt-i G G{X^^] 1) 
form a single chain, i.e., they satisfy the following properties: 

i) For each j ~ 1, 2, .... t — 2, (tij, Wj+i) e E. 

ii) For all 1 < j, fc < t - 1, {vj.Vk) i E for \j - A:| > 1, 

where E denotes the set of links in G{xi'-^^ \ 1). Then the t-th 
order cluster coefficient C^'^^ is defined to be the conditional 
probability that a node vt is adjacent to at least one of the 



nodes U2,...,wt-i, given that vt is adjacent to vi (averaging 
over all the possible positions ^.y^, . . . ,Xt,j ^ in R'^ of the 
points V2, ■ ■ ■ , wt-i satisfying conditions (i) and (ii)). 

According to the above definition, Cg''^ — C^'^K To calculate 
Cf'^^ is difficult in general. However, Cg''"', the cluster coeffi- 
cient for d-dimensional Poisson random geometric graphs can 
be computed as [15] 



(,) _ 3 r(^) 

' ^/^^(^) 



Tr/3 



l'' 0d0. 



(7) 



Among higher dimensional random geometric graphs, the 
three-dimensional case is of practical interest (as for sensor 
networks in the deep sea). Using the duplicate formula for the 
Gamma function, we can derive from (Ell as [15] 



(3) 



3/2 



2 ^/^.^,^(^ + i)V4 



0.4688. 



(8) 



Although it is difficult to obtain closed-from expressions 
for C['^\ t > 4, we are able to compute them by numerical 

(2) (2) 

integration. For example, we obtain CI ' « 0.6012 and Cf' « 
0.6179. We also note that < C^'^'' < 1, for all t > 3, since 
Cj"*^ is a (nonzero) conditional probability. 

C. Critical Density for Random Geometric Graphs 

Let H^^Q — 1-6^^ U {0}, i.e., the union of the origin and 
the infinite homogeneous Poisson point spatial process with 
density A in R''. Note that in a random geometric graph 
induced by a homogeneous Poisson point process, the choice 
of the origin can be arbitrary. 

Definition 3: For G{Ti!'^\^]l), the percolation probability 
Poo (A) is the probability that the component containing the 
origin has an infinite number of nodes of the graph. 

Definition 4: For G{li}^\^] 1), the critical density (contin- 
uum percolation threshold) Ac''' is defined as Ac'''' = inf{A > 
0:poo(A) >0}. 

It is known from continuum percolation theory that if A > 
X^f' , then there exists a unique connected component con- 
taining 0(n) nodes in G{xli'^] 1) a.a.s^This large connected 
component is called the giant component [7]. 

III. New Lower Bounds on the Critical Density 

A fundamental result of continuum percolation states that 
< X^f"^ < oo for all d> 2. Exact values for Ai'*' and poo(A) 
are not yet known. For d ~ 2, simulation studies [11] show that 

(2) 

Ac ~ 1.44, while the best analytical bounds obtained thus 
far are 0.696 < A^^' < 3.372 [6], [7]. Recently, in [12], the 

(2) 

authors reduce the problem of characterizing Ac to evaluating 
numerical integrals, and they obtain numerical bounds 1.435 < 

^We say f{n) = 0{g{n)) if there exists no > and constant cq such that 
fin) < co5(n) Vn > hq. We say /(n) = 0(5(n)) if g(n) = 0(f(n)). 
Finally, we say f(n) = B{g(n)) if /(n) = 0{g(n)) and f(n) = n{g(n)). 



Xi^^ < 1.437 with confidence 99.99%. Unfortunately, these 
bounds are not in closed form and are restricted to the two 
dimensional case. In the following, we present an analysis 
which combines a technique used in [13] (for random graphs) 
and the clustering effect in random geometric graphs to obtain 
a new mathematical characterization of the critical density Ac'''' 
for d> 2. This analysis yields a new class of improved lower 
bounds for A^''' . In particular, they yield the tightest analytical 
lower bounds known to date. 

Theorem 1: Let /i be the mean degree of G(A'i''';l), where 
d> 2. For any given integer t > 3, if 



1 



l~C 



where C^''' is defined by Definition 2, then the largest com- 
ponent of G{xi'^''; 1) has at most a Inn nodes a.a.s., where a 
is a positive constant. 

Before giving the proof, we define the diameter of a node 
with respect to a subgraph. 

Definition 5: Given a graph G — {V, E), for any subgraph 
C — {V',E') C G and a node u e V, define the diameter 
of node u with respect to G' — {V , E') as 



diam(M, G') = max{(i(u, w)}, 



(9) 



where d{u, v) is the distance between u and v, measured by 
the length of the shortest path between u and v in terms of 
the number of links. 

Note that the diameter of graph G is the maximum of 
the diameters of all nodes with respect to graph G, i.e., 
diam(G) = max„gG{diam(M, G)}. Another useful fact for 
the following proof is that for a random geometric graph 
G(A'„; 1), if diam(w, G') < c, then the Euclidean distance 
between u and any node v in G' — {V ,E') is no more than 
c, i.e., ||X„ - xJl < c, for all vC,V' 

Proof of Theorem |7] Let /i — ^^\d) , where < e < 1. 
For simplicity, let p denote the probability that there is an 
link between two nodes, i.e. p — Punk given by (O, so that 

[n - l)p = ^. 

We consider an arbitrary node (with fixed label and random 
position) V € G(A'i''\ 1) and study the following "active- 
saturated" process. For i = 0, 1, 2, . . ., let Ai denote the set 
of "active" nodes, and Si denote the set of "saturated" nodes, 
starting with — {v}, So ~ 0. At (i+l)-th step, we select an 
arbitrary node u from Ai and update the active and saturated 
sets as follows: 

A,+i = {A,\u) U {N, n (A, U S,Y), S,+i = 5*, U u, 

where Ni is the set of neighbors of u. In other words, at each 
step we move a node u from the active set to the saturated set, 
and at the same time move to the active set all the neighbors 
of u which do not currently belong to the active or saturated 



nodes Wj_^^l, 



set. In this manner, we can go through all the nodes in u's 
component, represented by Ty, until Ai = 0. 

Let Fi+i be the number of nodes added to Ai at step i + 1: C^"^^ (n, Uj-^-l,Wj_^_l 



in Sj that forms a single chain with 
node Wj+i (i.e., satisfies condition (i)-(ii) of Definition 2). Let 

be the conditional probability 



3 + 1' 



Y, 



i+l 



\N, n {A, U 5,)' 



(10) 



Note that \Si\ = i for i < |r„|. 

We say a sample graph G„ of G{X^'^ ; 1) is good, if there is 
no component having size strictly larger than In n, or for 
any node v e G„ with \Ty \ > Inn and any sequence of 
the "active-saturated" steps starting at v, there exists a bounded 
k' such that 

Vj > k' and Vu G Aj, diam(u, 5'^ U it) > t - 2, (11) 

Ck' = \Ak' U Sk'\ < {\iin)'^\ (12) 

Note that fc' depends on n, the sample graph G„, the node u e 
G„ and the sequence of the "active-saturated" steps starting at 
V. Let %i be the collection of all good sample graphs with n 
nodes. We will show later that with probability 1, there exists 
a uniform bound /cq < oo such that k' < ko for any n, any 
G„ e Tn and any w e G„ € T^j with |r„| > ^^\nn and 
any sequence of the "active-saturated" steps starting at v. We 
assume that this holds for the moment. 

Now given G(A',1''''; 1) G 7^, consider an arbitrary node 
(with fixed label) v e G{X,['^^ ;!), if |r„| > ^Inn, the 



"active-saturated" process can sustain at least — 
k = ^Inn - fc' > 0. 



In n steps. 



Because \Si\ — i, we have 

|r„| > A: + fc' ^ \Ak+k' U Sk+k' \>k + k'. (13) 

Since G{Xn'^'^;l) e 7^, with probability 1, there exist con- 
stants fco, and fc' < fco satisfying condition (fTTT i and (fTSI i. By 
the definition of c^', ( flj] ) is equivalent to 

k+k' 

|r,|>fc + fc-'^ J2 Y^>k + k'-ck'. (14) 



=l+fc' 



Therefore, 



Pr{|r^| > fc + fc'} 




> fc + fc' 



Ck' 



(15) 



We now bound the RHS probability. Consider the "active- 
saturated" process after fc' steps. For all j > fc', we have 
diam(M, Sj U u) > t — 2, Vii G Aj. At each step i, we move 
Yi nodes to the active set. Suppose at the {j + l)-th step, 
J2i=i — ™' there are n — 1 — m nodes remaining. 
Now suppose we move node Uj+i from the active set to the 
saturated set, then Y,+i is the number of nodes adjacent to 
Uj+i but not in Aj or Sj. Since all nodes in Aj are adjacent 
to some node in Sj, Yj+i is also the number of nodes adjacent 
to Uj+i, not in Sj and not adjacent to any node in Sj. Since 
diam(Mj+i,5j U Uj+i) > t ~ 2, there exists a sequence of 

^We ignore integer constraints for convenience. 



that one of the remaining n~l~m nodes, w, is adjacent to at 
least one of the nodes w]+i^ w*j+i given that w is adjacent 
to Uj+i. Then the probability that a node is adjacent to Uj+i 
and not adjacent to any w*^^ j * = 1, 2, f — 2, is (7„(uj+i) = 

where p„(uj+i) 



Pn ("j + l ) (1 - (n, Uj+i , wj+i , 

is the average probability that there is a link between node 
Wj+i and any other node. Since diam(wj+i,5j U Wj+i) may 
be larger than t — 2, and there are other geometric constraints 
for each of the n—1 — m remaining nodes, the probability of 
any one of the remaining n — m—1 nodes is adjacent to Uj+i, 
not in Sj and not adjacent to any node in Sj is less than or 
equal to q„(uj+i). 

Now Yj+i = ™ where Bj = 1 if node i is 

adjacent to Uj+i, not in Sj and not adjacent to any node 
in Sj, and Bi — otherwise. Note that these Bi's are not 
independent. Nevertheless, by the argument above, we have 

Pl{Bi + i = l\{Bi,B2,...,Bi) = {h,b2,...,bi)} < qn{Uj + l) 

for any (fei, &2, 6;) S {0; 1}' ™d ^ ^ n ~ 2 — m. For 
i = 1, . . . ,n — 1 — m, let B^ be i.i.d. Bernoulli(g„(uj+i)) 
random variables. By Proposition 1 in Appendix I, = 
™ is stochastically upper bounded by Y"_^-^ = 
^ Binom(n — 1 — m,q„(uj+i)), which is 
further stochastically upper bounded by ij+i ^ Binom(n — 
1, qn{uj+i)). Therefore, conditional on X]i=i ^ ™' for ^i^Y 
m < n—1, Yj+i is stochastically upper bounded by a random 
variable K,'+i with distribution Binom(n — l,q„(uj+i)). 

Using the same argument, we see that '}2!i=i+k' 
is stochastically upper bounded by X]i=^i+fe' ^/ = 
Zk ^ Binom(fc(n - l),g„(u„j)), where qn{um) = 

SUPl + fe'<j<fc+fc' (ln{u^). Thus, 

k+k' 1 

Pr<j Y,>k + k' -Ck'\<Vr{Zk>k + k' -Ck'}. 

i=l+k' \ 

(16) 

Let ^n{um) = {n — l)pn(um) be the mean degree of node 
Um- Since E[Zk] = k{n - l)qn{um.) = kfj,nium)il - 
Cj.'^\n,um,w'^, ...,wli)), 

Pr{Zk >k + k' -Ck'} 
= Pi{Zk > E[Zk] + k + k' - Ck' - E[Zk]} 
^ Pr{Zfc > E\Zk\ + fc<5„(M,„) + fc' - Ck'}. 

where 5„(um) = 1 - /Lt„(w„i)(l - cf\n,Um,w\^, w^)). 
Note that conditioned on G(A'i'''';l) G T^i, the node distri- 
bution may not be uniform. Nevertheless, we will show that 
Pr{G(A'i''^ 1) G T„} ^ 1 as n ^ oo. Hence the node 
distribution is uniform asymptotically, and p„(uj+i) -p. 



■' ""1+1 



firi{uj+i) n and G^^ •'(n, m^+i, wj^;^ 
n ^ oo. Thus 5n{um) — > e as n — > cx). By (l22T i. there exists 
< no < oo, such that for n > no, \Sn{um) — e| < f and 

kSn{Um) + fc' - Cfc' 



3 — 2e 

= Sn{u„i)—-^ Inn + {I - 6n{Um))k' ~ Ck' 

> in 71 — In n 3 

2e ^ ' 

> 0. 

By the Chernoff bound [13], for S > 0, 
Pi{Z > E[Z] +5}< exp 



2E[Z] + 2(5/3 



(17) 



Thus, for n sufficiently large 

Pr{^fc 5 

< exp < — 



Pr{Zfc > fc + fc' -Cfc/} 



k' 



2fc(l )) + 2(fc5„(i(™) + fc'-CfeO/3 



+l-Pr{G(A'W;l)^r„} 
= 0(n-^)+Pr{G(A'W;l)^T„}. 

To complete the proof, we show two facts: 

(i) Pr{G(A'i'^^ 1) e r„} ^ 1 as n ^ oo. 

(ii) With probability 1, there exists < oo such that 
k' < fco for any n, any G(A'^'''';1) e Tn and any 
V G G(A'^''''; 1) e 7^1, and any realization of the "active- 
saturated" process starting at v. 

To show (i), note that in G(X^^ ; 1), if there is no compo- 
nent having size strictly larger than a Inn, then it is good; 
otherwise, we prove that for any node v E 



exp 



k'^S-aiUmf + 2fc5„(ll™)(fc' -Cy) + {k' - CfcO 



fc(3 - 2S„(um)) + k' - Cfc/ 

Since k = ^^^Inrt — k', k' < kg and Ck' < (Inn) 3, as 
n oo, the RHS of ( fTSl ) has the same order as 

' 3r(3-2e)(5„(w 



with |r„| > 



In 71 and any realization of the "active- 



exp ■ 



< exp 



exp 



Inn 

2(5„(Um) SniUrn) 



e2(3-2,5„(u™)) 



3 - 2Sn{u,n) 
3^n(Mm)(2 - SnjUjn)) 
2(3 - 2(5„(l7„)) 

3(3-2e)5„(M™)2 

Inn - 



1 



Cfc' 



)]} 



35„(u„ 



Cfc-' Wl9) 



2e2(3-2(5„(u™)) ' 3-2(5„(m™) 
Now choose 7 > such that ''V^o^'J'^^^'^^} > i- Because 

' [3 — 2(c — 7)] 6 

Snium) e as 71 — > 00, there exists < ni < 00, such 
that for 71 > 7ii, \6n{um) — e| < 7- Then using e — 7 < 
5n{um) < e + 7, we can bound fT% by 



c exp 



3(3-26)(6-7)^ 
2e2[3-2(e-7)] 



Inn 



3(e + 7) 



3-2(e + 7) 



< c' exp ^ — - In n 



3(e + 7) 



rCfe' 



0(n-3), 



where 



exp 



3-2(e + 7) 
3(6-7)[2-(e + 7)] 



2[3-2(6-7)] 
By O-fEOl), for any ai'biti-ary node v £ G(A'i''^ 1) e T„, 

Pr{|r,| > Inn} =0(n-3). (21) 



saturated" process starting at v, there exists a bounded fc' such 
(1^)' that (fTTT l holds a.a.s. Suppose for some w e G(A',1'''; 1) with 
\^v\ > ^^^^Inn and a realization of the "active-saturated" 
process starting at v such that for any k < |r„| — 1, there 
exists a step j,k < j < |r„| — 1, and w E Aj, such 
that diam(w,S'j U w) < t - 2. Since |r„| > ^^^Inn, as 
n — > 00, |r„| 00, the "active-saturated" process can go on 
forever. Thus, Sj asymptotically contains an infinite number 
of nodes. Since diam(it;, Sj U w) < t ~ 2, all the nodes of 
Sj lie in a ball centered at w with radius t — 2 (by the 
argument immediately following Definition 5), which occurs 
with probability approaching as 71 ^ cx), by Lemma 1. 
Now suppose k' satisfying (fTTT l exists but is unbounded as 
71 ^ CX3. Then let k' be the smallest k' satisfying (fTTT l. For 
step j = k' — 1, there exists at least one node w E Aj such 
that diam(7i;, Sj Uw) < t — 2. Then, arguing as before, we can 
show this happens with probability approaching as n ^ cx). 

Next we show that Ck' < (Inn) 3 a.a.s. We know Ck' — 
\Ak' U Sk'\ = J2i=i^i- Using an argument similar to that 
leading to ( fTST l. we see that the RHS is stochastically upper 
bounded by a random variable with distribution Binom(fc'(n — 
Assuming fact (ii), which is shown below, the RHS is 
further stochastically upper bounded by a random variable Ck' 
with distribution Binom(fcon,p). Applying the Chernoff bound 
(fTTT l. we have, for sufficiently large n, 

Pr{cfc' > (ln7i)3} <Pr{cfc' > (Inn)^} 

= Pr{cfc' > ko^ + (Inn) 3 - kofi} 

[(In n)^ — /co/i]^ 1 



Cfc' 



(20) 




2fco/i + I [(In n) 3 - kofi] 



Set a = The probability that random geometric graph 

G(A'^'^'';1) has at least one component whose size is no 
smaller than a In n is 

Pr{3t; E G{X^'^^; 1) : |r„| > a In 71} 
= Pr{3t; E G(A'^'^); 1) : |r,| > aln7i|G(4'^); 1) E T„} 
•Pr{G(A-('^);l)e7;,} 

+Pr{37; E G{Xi''^; 1) : |r,| > alnn|G(A'W; 1) i %,} 
•Pr{G(A'W;l)^r„} 
<nPr{|r,| > alnn|G(4'^); 1) E r„}Pr{G(4''); 1) E %,} 



a.a.s. 



(22) 



we show fact (ii). Suppose that for each 
., there exists ii^ and v E G{X^^]\) E 7^. 



Finally, 

I = 1,2, 

with Ui > \T^\ > 2^^1nni and a realization of the "active- 
saturated" process starting at v such that \Ty\ ~ 1 > k' > i. 
As 1 ^ 00, fc' ^ 00, \Ty\ — + 00, and rii 00. However, this 
holds with probability 0. This completes our proof. □ 



Note that by either ignoring border effects or by taking 
n 00, C^'^^ (defined for G{X^'^^; 1)) is also the t-th order 
cluster coefficient for the infinite Poisson random geometric 



graph G{T-L^^^ ; 1). Thus by Theorem[T] we obtain the foUowins 
important corollaries giving new improved lower bounds on 
the critical density for d-dimensional Poisson random geomet- 
ric graphs. 

Corollary 1: Let /ii*^^ and A^'''' be the critical mean degree 
and critical density for G{Ti^^''; 1), respectively, where d> 2. 
Then for all t > 3, 

Mi'^> and A(^)>— —. (23) 



(d) 



I -a 



id) 



Proof : Follows immediately from Theorem [T] and □ 



In particular, for two-dimensional Poisson random geomet- 



(2) 



ric graphs, substituting V^'^^ = tt and C; 
0.5865... into ( |23] ), we have the following corollary. 

(2) 



3^3 



Corollary 2: The critical mean degree /ic and the critical 



density A, 
0.7698 



(2) 



forG(7^f^l) satisfy fi, 



(2) 



> 2.419. 



and A^^^ > 



(2) 

Note that if we use high-order cluster coefficients ,t > 



(2) 

4, computed by numerical methods, e.g., C4 w 0.6012 and 



C. 



(2) 



0.6179, we can obtain further improved (approximate) 



lower bounds: > 2.617, and Ac > 0.8830 

By applying Cf = 0.4688... given by ([8]), we have 

Corollary 3: The critical mean degree /ii^^ and the critical 



density A, 
0.4494.... 



(3) 



for G(Hf ^ 1) satisfy > 1.412... and A^^^ > 



C3) 



.(3) 



This lower bound is close to the known results obtained 
by simulation — 0.65 [16], and it is the best known analytical 



lower bound on A, 



(3) 



IV. Conclusion 



We have established a new class of analytical lower bounds 
on the critical density xlf'' for percolation in d-dimensional 
Poisson random geometric graphs. These analytical lower 
bounds are the tightest known to date, and reveal a deep 
underlying relationship between the cluster coefficient and the 
critical density in continuum percolation. 
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Appendix I 

Proposition 1: Suppose random variables Xi,i = l,...,m 
satisfy the following conditions: (i) Pr{Xi > x} < Pr{X^ > 
x}, Mx- (ii) Vy{Xi > x\Xi = xi,...,Xi^i = xi^i} < 
Pr{X^ > x}, yx.xi, ...,xi^i\ (iii) X^ ,i — l,...,m are 
independent of each other and of X^, i = 1, to. Then 



Pr 5]X,>z <Pr^X+>^ 



Vz 



Proof: It suffices to show the result for m = 2. Since 
Pr{X2 > y\Xi - xi} < Pt{X+ > y\Xi = x,}, Va;i,y, 
Pt{Xi + X2 > z\Xi = a;i} < Pr{Xi + X+ > z\Xi = 
xi}, Vxi,z, and thus 

Pr{Xi +X2>z} 

= E.iPr{^i=2;i}Pr{Xi+X2>z|Xi-xi} 

< = xi}Pr{Xi + X+ > z\Xi = xi} 
= Pi{Xi + X+ > z} 

= E4 Pr{^2+ = }Pr{^i + 4 > ^1^2+ - 4} 
= E4 Pr{^2+ = ^^}Pr{^i >^-4} 

< Pi-{^2+ = 4 }Pi-{^i+ >^~4} 

= E4 Pi-{^2+ = }Pr{^i+ + X+ > z\X+ = x+} 
= Pi{X++X+>z}. 

□ 



